#|# output-location: slide # Puts output on the next slide
results_pm_copd <-
healthiar::attribute_health(
erf_shape = "log_linear",
rr_central = exdat_pm_copd$relative_risk,
rr_lower = exdat_pm_copd$relative_risk_lower,
rr_upper = exdat_pm_copd$relative_risk_upper,
rr_increment = 10,
exp_central = exdat_pm_copd$mean_concentration,
cutoff_central = exdat_pm_copd$cut_off_value,
bhd_central = exdat_pm_copd$incidents_per_100_000_per_year/1E5*
exdat_pm_copd$population_at_risk
)